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We study doping dependence of interlayer hoppings, t±, in multilayered cuprates with four 
or more CuC>2 planes in a unit cell. When the double occupancy is forbidden in the plane, an 
effective amplitude of t± in the Gutzwiller approximation is shown to be proportional to the 
square root of the product of doping rates in adjacent two planes, i.e., t c ^ oc t ±1/6182, where 
Si and 82 represent the doping rates of the two planes. More than three-layered cuprates have 
two kinds of Cu02 planes, i.e., inner- and outer planes (IP and OP), resulting in two different 
values of t e ±, i.e., tf±i c< i±V<5ip<5ip between IP's, and t e ± 2 oc i±V<5ip(5op between IP and OP. 
Fermi surfaces are calculated in the four-layered t-t'-t"-J model by the mean-field theory. The 
order parameters, the renormalization factor of t±, and the site-potential making the charge 
imbalance between IP and OP are self-consistently determined for several doping rates. We 
show the interlayer splitting of the Fermi surfaces, which may be observed in the angle resolved 
photoemission spectroscopy measurement. 
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1. Introduction 

Cuprate superconductors have layered structure of 
CuC>2 planes, which constitute conducting blocks 
separated by charge-reservoir blocks. Multilayered 
high-T c cuprates, e.g., Ba2Ca 3 Cu408(Oi_, / F, / )2 1-3 and 
HgBa2Ca4Cu50j,, 4 ~ 6 have two kinds of CuC>2 planes 
in one conducting block; the pyramidally-coordinated- 
outer-planes (OP's) and the square-coordinated-inner- 
planes (IP's). Due to spatial variation of Madelung po- 
tential, 7-10 the carrier density in the OP is generally dif- 
ferent from that in the IP. We call such an inhomoge- 
neous charge-distribution 'charge imbalance'. 

The charge imbalance has been observed by the nu- 
clear magnetic resonance (NMR) measurement 11-17 and 
the X-ray diffraction measurement with the bond- valence 
analysis. 18-20 As regards hole-doped systems, the NMR 
studies have found that a hole concentration in the OP, 
Sop, is larger than that in the IP, <5ip. A difference of 
doping rate, AS = Sop — Sip, increases with the number 
of Cu02 planes in one conducting block, n. 17 Although 
the five-layered cuprates have three IP's, the charge im- 
balance among them is small compared to AS. 

It is known that the superconducting critical temper- 
ature, T c , increases with n and saturate at n = 3. 18 ' 19 
Thereafter, T c decreases with n, while the charge imbal- 
ance increases. Considered such an n-dependence of T c , 
the charge imbalance seems to suppress T c for n >3. Si- 
multaneously, the charge imbalance induces some inter- 
esting phases, e.g., two kinds of superconducting (SC) 
gaps 16 and coexistence of the SC- and antiferromagnetic 
(AF) states. 21,22 In the former case, the OP's have a 
large SC gap compared to the IP's. In the latter case, 
the SC OP's are separated by the AF IP's, and a long- 
ranged Josephson coupling through the AF planes stabi- 
lizes the SC state as a bulk. 23 The suppression of T c may 
be partly caused by these inhomogeneous phases. 



On the other hand, some theoretical studies have been 
presented that an interlayer hopping plays an important 
role about the n-dependence of T c . 23-26 The interlayer 
hopping is indispensable for the bulk SC state, even if 
each Cu02 plane shows superconductivity. 27-33 Here, we 
note two kinds of interlayer hoppings in the multilayered 
cuprates. One is a coherent hopping between adjacent 
planes within the conducting block, and the other one 
is an incoherent hopping through the charge reservoir 
block. The 'coherent' means the conservation of momen- 
tum. Although both contribute to stabilize the bulk su- 
perconductivity, it is the former interlayer hopping that 
gives rise to the n-dependence of T c . Below, we will fo- 
cus our attention on the coherent hopping within the 
conducting block and simply call it 'interlayer hopping'. 

An amplitude of the interlayer hopping, t±, has been 
observed in the bilayer cuprates, e.g., Bi2Ba2CaCu20j / , 
by the angle-resolved-photoemission-spectroscopy 
(ARPES) measurement. 34 ' 35 The experimental estima- 
tion of t± in the bilayer one is about 50 meV, 34 which is 
much smaller than 150 meV obtained by the bilayer LDA 
calculations. 36 On the other hand, the bilayer Hubbard 
model predicted a similar amplitude of 40 meV. 37 This 
is because a strong on-site Coulomb repulsion in the 
Cu02 plane substantially reduces i^. 38 ' 39 Furthermore, 
the interlayer hopping depends on the doping rate as 
well. For example, in the bilayer case, the splitting of 
Fermi surface (FS) is clear in the overdoped region, 
while it is not observed in the underdoped region. 34 ' 35 
Therefore, as regards the multilayered cuprates for 
n >4, the relation between the interlayer hoppings and 
the Fermi surface splittings are not so obvious. One 
can suppose that the charge imbalance gives spatially 
inhomogeneous t±'s, i.e., t± between IP's is different 
from that between IP and OP. Such inhomogeneous 
interlayer hoppings lead to unexpected FS splittings, 
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which will tell another key factor of n-dependence of 
T c for n > 3. 

In this paper, wc study doping dependence of t± to 
show the FS splitting in the multilayered cuprates. The 
effective t± renormalized by the on-site Coulomb re- 
pulsion is calculated in the Gutzwiller approximation 
(GWA). 40-53 It has been shown that the GWA gives a re- 
liable estimate of the variational energy for the metallic- 
and SC states, 48-50 while extensions of the GWA are 
necessary for the AF state. 51-53 Below, we consider the 
metallic CuC>2 planes without the AF order. When the 
double occupancy is forbidden in each plane, we have 
found that the renormalization factor, is proportional 
to the square root of the product of the doping rates in 
two planes, resulting in two different values of q±, i.e., 
3_L,i oc V^ip^ip between IP's, and q±^ oc V^ip^op be- 
tween IP and OP. Taking account of q± t \ and q± t 2, FS's 
are calculated in the multilayered t-t'-t"-J model in the 
mean-field theory. The order parameters, the renormal- 
ization factors, and the site-potential making the charge 
imbalance are self-consistently determined for several 
doping rates. We find that the FS splittings are not al- 
ways well separated, e.g., only two FS in the four- layered 
system. 

The rest of this paper is organized as follows. In §2, 
wc briefly summarize the Gutzwiller approximation for 
the interlayer hopping between two planes with different 
doping rates. In §3, the FS's are calculated for the four- 
layered case of the multi-layered t-t'-t"- J model including 
the inter-layer hopping obtained in the previous section. 
Summary and discussions are given in §4. 

2. Gutzwiller Factor of Interlayer Hopping 

We consider two planes with the charge imbalance 
given by 



H — Hn + H\y + H±, 
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where cf^ {cfl) l& ^ ne electron creation (annihilation) 
operator with spin a at the i-th site in the l-th plane, and 

n\2 = c i^ c( ia and = n i\ + n ii- Hopping integral in 
the l-th. plane and a strength of Coulomb repulsion are 
denoted by and U, respectively. Several cases of the 
charge imbalance are schematically shown in Fig. 1. The 
hole density is denoted by 5, and the density of double 
occupancy is indicated by d. 

The Gutzwiller's wavefunction is assumed for the 
ground state as 



(5) 



where P n is the projection operator onto the subspace of 
fixed densities, and n^ 2 \ The Gutzwiller's projection 



operator, P g , with a variational parameter, gi, is given 

by 



(6) 
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and \ipo) is a ground state wavefunction for U = 0. 
The ground state energy, E g , is evaluated as 



(7) 



In the GWA, P g is replaced by classical weighting fac- 
tors, which is called the Gutzwiller factor. Accordingly, 
Eg is calculated by the expectation value of an effective 
Hamiltonian as, 



E s = (Vol-ffefflV-o) +U(di +d 2 ), 
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where di is a density of doubly occupied sites in the l-th 
plane. The Gutzwiller factors of the intraplane- and the 
interlayer hoppings are denoted by q^ and q±, respec- 
tively. 

In the original GWA, only the site-diagonal expecta- 
tion values, (n-o-), are considered. In contrast, the inter- 
site correlation, {c\ a Cj a ), becomes important to obtain 
qf^ in the AF state 44, 51-53 having the two sublattices. 
Then, it is interesting whether the intersite correlation 
for q± is important in the paramagnetic state with charge 
imbalance. In this case, however, the intersite correlation 
is negligible for q±. Details are summarized in the Ap- 
pendix. 
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In the way of original GWA, q^ is given by 
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Even in the site-diagonal approximation, q± is not so 
obvious, since the two planes have different carrier den- 
sities. 46 ' 47 The renormalization of the interlayer hopping 
from the Z-th plane to the m-th one is different from that 
of the conjugate process as, 
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The GW factor for the interlayer hopping can be ob- 
tained as, 
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If one considers two hole-doped planes without double 
occupancy, i.e., d\= d 2 = 0, the interlayer GW factor is 
given by 

91 - ( ( i + <.)(i + w j • (17 » 

where (5; = 1 — 2ni. 

If the first plane is a hole-doped one without double 
occupancy and the 2nd plane is an electron-doped one 
without hole, i.e., 1 — 2n 2 + d 2 =0, the interlayer GW 
factor is given by 

^ = ( „.t hd ,\,X ", us, 



(i + £i)(i + d 2 )y ' 

where 1 — n 2 = n 2 — d 2 . 

3. Fermi Surfaces in Four-Layered Systems 

3.1 Multilayered t-t'-t"-J model 

To examine the FS's in four-layered cuprates, we adopt 
multilayered t-t'-t"-J model 54 ' 55 composed of two IP's 
and two OP's as shown in Fig. 2. The Hamiltonian is 
given by 

H = H\\+H ± +H W , (19) 
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.) is the annihilation opcr- 
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ator of an electron in the l-th. layer with spin a at 
site i with the constraint of no double occupancy, and 



,(0 



hi 



and £p arc the charge and the spin 



operators, respectively. 

In eq. (20), the summations (ij)ut, (ij)2nd and {ij) 3r d 
run over first, second and third nearest-neighbor sites, 
respectively. The values of the parameters are as fol- 
lows; J = 0.14 cV, t/J = 2.5, t'/J = -0.85 and 
t"/J = 0.575. 54 

The inter-layer hopping in eq. (21) has the dispersion 
relation, 30 

e±(k) = ~ {cos{k x ) - cos(k y )) 2 , (24) 

where t± is the amplitude without renormalization and 
is set to t±/J= l.O. 30 - 34 - 37 

The chemical potential, fi, is introduced to fix the total 
density as, N — N\p + Nqp- The charge imbalance is 
measured by a difference of density in the IP and that 
in the OP, i.e., SN = N\p — Nop. SN is determined by a 
site-potential in eq. (22), W, and the effective interlayer 
hopping, t e g = q±t±. In this paper, a value of W is self- 
consistently determined by fixing SN . 

The order parameter given by 



is introduced to the fourth term in eq. (20). The bond 
directions are indicated by r = x, y. The parameters, /x, 
W, and Xr ) are self-consistently determined by numer- 
ical calculations. 

We adopt the Gutzwiller approximation shown in §2 
to take the constraint of no double occupancy in the 
Hamiltonian (19). As a result, the electron creation- and 
annihilation operators arc released from the constraint, 
and the parameters are changed as follows: 

t -+ qfh, t'^q^t', f'^qfh", (26) 

J -> qfj, (27) 

tj_ — > tj_i = q±it_L, t^ 2 = q± 2 tj_, (28) 

with 
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If the IP becomes electron-doped, Sip is replaced by dip. 
The effective amplitudes of interlayer hoppings are de- 
noted by tj_i = (7_Litj_ between IP and OP, and t± 2 = 
q± 2 t± between IP's. 



(0 
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3.2 Fermi surfaces 

Fermi surfaces (FS's) of four-layered cuprates are 
shown in Fig. 3 for (a) # av =0.15, 5ip=0.15, <5op=0-15, (b) 
(5 av =0.15, 5ip=0.11, (5 OP =0.19, (c) <5 av =0.05, <5 IP =0.01, 
(5 OP =0.09, and (d) <5 av =0.00, <5 IP =-0.10, (5 OP =0.10. The 
averaged hole density, <5 av , is defined as S av = (Sjp + 

<y OP )/2. 

In the case of no charge imbalance given by (5ip=0.15 
and <5op=0.15, the FS does not split along the nodal di- 
rection as shown in Fig. 3 (a), since e±(k) in eq. (24) 
becomes zero in this direction. If one assumed a 'fc- 
independent' interlayer- hopping the FS could split 
along the nodal direction as well. However, by taking 
account of the experimental results about the bilaycr 
cuprates, an amplitude of tj_ must be small compared 
tot ± . 34 - 35 

On the other hand, in Fig. 3 (b), the FS along 
the nodal direction splits due to W that induces the 
charge imbalance given by <5ip=0.11 and <5op=0.19. The 
Fermi points in this direction are determined by — 
ej(fe) ± W/2 - fj,, where ej(fe) = -4i; cos/c - 4iJ(cos k) 2 - 
At['cos2k- (3 / '2) Jixi cos k, U = q[ l) t, t[ = qfh', t'{ = 
qf't", and J; = q*pJ. The width of the splitting is 
roughly estimated as J2i W/{8t t + 3 J/Xi) for t' = t" = 0, 
fj, = and W/ti,W/(JiXi) <C 1. Without the charge im- 
balance, the FS does not split along the nodal direction, 
since W — > and e±(k) in eq. (24) becomes zero in this 
direction. Therefore, one can say that the FS splitting 
along the nodal direction stands for the charge imbal- 
ance. 

In both cases, Figs. 3 (a) and (b), the systems have 
enough hole density as <5 av =0.15, which splits the FS 
into four. Two of them look like a hole-like FS enclos- 
ing (ir, 7r) point, while the other two enclose (0, 0) point 
like an electron-like case. Such two different types of FS 
originate from a fact that the dispersion relation is flat 
around the antinodal points, i.e., (±7r, 0) and (0, ±tt). 
As a result, a tiny increase of electron density leads to 
a large shift of FS. Particularly, since an amplitude of 
e±(k) is maximum at the antinodal points, the interlayer 
splitting is large and then results in the electron- and 
the hole-like FS's. Note that there is no electron-doped 
Cu02 plane in Figs. 3 (a) and (b). 

By decreasing the averaged hole density as <5 av =0.05 
and imposing the charge imbalance given by <5ip=0.01 
and <5op=0.09, the number of the FS's seems to ap- 
proach to two as shown in Fig. 3 (c); one is hole- like 
FS and the other one is electron-like FS. The split- 
ting of the two hole-like (electron-like) FS's is deter- 
mined by ij_i/JipXip, which is proportional to Sip, and 
i_i_i/t_i_2 = V ^ip/^op- 

So far, we have discussed the hole-doped cases given 
by <5 av > 0. Assumed that all apical oxygens are ideally 
substituted by fluorine like Ba2Ca3Cu40sF2, the nomi- 
nal valence of Cu becomes +2. In such a case, we can ex- 
pect the half-filled system as <5 av =0.0, but the system has 
the charge imbalance given by S\p=-dip and <5op = dip- 
Note that the IP becomes electron doped. In Fig. 3 (d), 
we show FS's in such a case with Sop=dip = 0.10. The 
FS's are split into four as is the case of Fig. 3 (b). 



The splitting or the degeneracy of FS's is explained by 
a simple energy diagram shown in Fig. 4. Without the 
interlayer hoppings, the degenerated OP's are located at 
W in energy above the degenerated IP's. By considering 
the inversion symmetry, each pair of planes is classified 
into even- and odd wavefunctions given by OP cvcn .odd = 
(OPi±OP 2 )/V2 and IP C vc„,odd = (IPi±IP 2 )/V2. The 
IP cvcn and IPodd are separated by t±i, which finally 
leads to the splitting between the hole-like (electron-like) 
FS's. The remaining degeneracy is lifted by t±2- Conse- 
quently, we obtain four separated FS's, |1), |2), |3), and 
|4), from (0, 0) to (it, tt) points as indicated in Fig. 3 
(b). Each wavefunction is given by the linear combination 

as, |1) EE a o 0P od d+/3oIPodd, |2) = a c OP even +/3cIPcvcn, 

|3) ee <OP odd +#,IP odd , and |4) ee <OP ovcn +/3aP cvcn . 
The splitting between the hole-like (electron-like) FS's in 
Figs. 3 (a)-(c) is identified with that between |3) and |4) 
(|1) and |2)). Although the hole-like FS's are mainly com- 
posed of IP's wavefunction, the OP's wavefunction also 
contribute to the hole-like FS's. As a result, for example, 
a ratio of the area of FS |4) to the Brillouin zone is equiv- 
alent to the electron density given by a' 2 nop+/3' 2 nip, 
where nop and n\p are the electron densities in the OP 
and the IP, respectively. 

4. Summary and Discussion 

In this paper, we have examined the charge imbalance 
effects on Fermi surface (FS) splittings of multilayered 
cuprates by the multilayered t-t'-t"- J model in the mean- 
field theory. Amplitude of interlayer hoppings is renor- 
malized by the electron-electron interaction in the Cu02 
plane. In order to take such a correlation effect in the cal- 
culation, we have extended the Gutzwillcr approximation 
to the interlayer hopping and have shown its Gutzwillcr 
factor. 

It is found that the splitting of FS's along the nodal 
direction stands for the charge imbalance, which leads to 
the inhomogeneous interlayer hoppings. Considered no 
double occupancy in the Cu0 2 plane in the four-layered 
system, the Gutzwiller factor, q±, is proportional to the 
square root of the product of the doping rates in neigh- 
boring two planes, resulting in two different values of 
q±, i.e., q±_i oc V^ip^ip between inner planes (IP's), and 
q±,2 oc ^/5~lp5~op between IP and outer plane (OP). With 
decreasing a hole density in the IP, four FS's in the four- 
layered system degenerate into two due to q± t i — > 0. 

In this study, we have not considered the antiferromag- 
nctic (AF) long- nor short-ranged order, which induces 
gapped behavior in the single-particle spectral function, 
A{k,uo), near the antinodal region. 56 ' 57 The gapped be- 
havior appears in the FS observed by ARPES as an 
electron pocket. 58 ' 59 Considered the charge imbalance 
in the underdoped region, the IP's may have the AF 
short- or long-ranged order, and are connected to the 
paramagnetic- or the superconducting (SC) OP's by the 
interlayer hoppings. To compare our results with exper- 
imentally determined FS's, it may be necessary to take 
the AF order in the calculations but leave them as future 
issues. 
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Appendix: Gutzwiller Approximation with 
Inter-Site Correlation 

The derivation of eq. (16) is summarized in this Ap- 
pendix. We apply the method developed by Ogawa and 
Kanda (ref. 44) to include the intersite correlation into 
Q.L- 

First, we consider two planes with L sites per plane. 
The number of up- and down spins in Z-th plane are de- 
noted by Mi and Ni, respectively. Given constants, L, 
Mi and Ni, satisfy the following relations as 

L = Di + Ai+Bi + E h (A-l) 

Mi = Ai+Di, (A-2) 

Ni = Bi + Di, (A-3) 

where Di and £/ are the number of the double-occupied 
and the vacant sites in the Z-th plane, respectively. Ai 
(Bi) means the number of single-occupied sites in the 
l-th plane with up (down) spin. Next, one divides the 
system into pairs of two sites, in which one belongs to 
1st plane and the other one does to 2nd plane. Each pair 
has various configurations, e.g. one is up-spin site and 
the other one is doubly occupied site. The number of 
pairs with Si-site in the 1st plane and S2-site in the 2nd 
one are denoted by y SlS2 (sj = Si/L,Si = D h A h B u Ei). 
These quantities, y slS2 , satisfy the following relations, 
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Pi/L and qi = Qi/L. The parameters, pi, qt, 



and y Sl s 2 , are determined by optimizing the wavefunction 
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- = (cffM^o = (c^c^) . (A26) 

The number of configurations that equally contribute to 
the summation in eq. (A- 14) is given by 



(A-27) 



In the limit of infinitely large system, the summation is 
replaced by dominant terms, and eq. (A-14) is approxi- 
mated as 
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where "(0)" indicate the dominant term. Below, we will 
abbreviate "(0)" to simplify the notation. By optimizing 
the following function, 
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+ 



2D X 2D 2 

9i 92 



LI 



llsi,S2 ^ s l s 2 ' 



detWdct V 



hiA^Y^-Sx), 

si=di,oi,6i,ei S2 

+ J2 mA S2 (^y siS2 -S 2 ), 

S2=d 2 ,a 2 ,b 2 ,e 2 Si 

+ E ln < >S2 (E y ^^-^) 
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+ S kfti^E^-G*)' (A-29) In the case of £7 > oo 

Sl ,s 2 ,i 

one can obtain the dominant terms in a factorized form 
as 



1/2 



Y Sl s 2 ^S\^Si^ L S\Sii^S\S'2- (A'30) 

The parameter A's introduced in eq. (A-29) as the La- 
grange multiplier are determined by solving the simulta- 
neous equation as, 

D\ = A dl A d2 aif3i + A dl A a2 ai(3 2 + A dl A b2 a 2 (3x 

+ A dl A e2 a 2 f3 2 , 



(A-31) 



A x = A ai A d2 ax/3 3 + A Ql A a2 ai/3 4 + A ai A b2 a 2 p 3 
+ A ai A e2 a 2 (3 4 , 

Bx = A bl A d2 a 3 f5x + A bl A a2 a 3 (3 2 + A bl A b2 a 4 (3x 
+ A bl A e2 a 4 (3 2 , 

Ex = A ei A d2 a 3 (3 3 + A ei A a2 a 3 [3x + A ei A b2 a 4 (3 3 



+ A ei A e2 a 4 (3 4 , 

D 2 = A dl A d2 ax(3x + A ai A d2 ax(i 3 + A bl A d2 a 3 (ix 
+ A ei A d2 a 3 (3 3 , 

A 2 = A dl A a2 ax(i 2 + A ai A a2 ax(3 4 + A bl A a2 a 3 f3 2 
+ A ei A a2 a 3 /? 4 , 

B 2 = A dl A b2 a 2 f3x + A ai A b2 a 2 (3 3 + A bl A b2 a 4 Px 
+ A ei A b2 a 4 f3 3 , 

E 2 = A dl A e2 a 2 p 2 + A ai A e2 a 2 p 4 + A bl A e2 a 4 f3 2 
+ A ei A e2 a 4 /3 4 . 
The optimizations about Di as dF/dDi = lead to 

(A-39) 



(A-32) 
(A-33) 
(A-34) 
(A-35) 
(A-36) 
(A-37) 
(A-38) 



2 _ A Dl A El 
.9; - 



Although there are several ways to get the Gutzwiller 
factor of hopping integral, we follow the way of Vollhardt 
(ref. 45). Thereby, the amplitude of interlayer hoppings 
are estimated as, 

y ai e 2 + 92y ai b 2 + giyd 1 e 2 + gxgiyd^ , , , (2)t (i), , x 



WMi(l — Wm 2 ) 



1(7 J<7 

(A-40) 



_ Veia 2 + giybia 2 + 92Vetd 2 + 929xyb 1 e 2 

wm 2 (1 — wmJ 



(A41) 

Finally, we find that the Gutzwiller factor of the inter- 
layer hopping is given by 

rm^cfM) (*\&»$\*) V /2 _ (A42) 



, , (A-43) 

We have numerically solved the simultaneous equation 
(A.31)-(A.38) as a function of u and v. It is found that 
the intersite correlations are negligible for small values of 
u and v about less than 0.2. In the multilayered t-t'-t"-J 
model introduced in §3, u and v are smaller than 0.1, and 
then the site-approximation can be justified. Equation 
(16) is obtained in the site approximation of eq. (A-43). 
In this case, u = v = and the following relations are 
obtained as, 



WtJ2)| 



12 
13 

14 

15 

16 
17 

18 
19 
20 
21 

22 

23 



9i 



di(l - 2ni + di) 



{ni 



A s , A s 



= L y SlS2 



SxS 2 
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Fig. 1. (a) equivalent two planes with hole-doping, (b) hole-doped 
two planes with charge imbalance, (c) the one is hole-doped and 
the other one is electron-doped. 
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Fig. 3. Fermi surfaces of four-layered systems with J = 0.14 
cV, t/J = 2.5, t'/J = -0.85, t"/J = 0.575, t ± /J = 1.0 
for (a) 5 av =0.15, 5 TP =0.15, 5 O p=0.15, (b) 5 av =0.15, 5 IP =0.11, 
5 O p=0.19, (c) 5 av =0.05, <5ip=0.01, 5 O p=0.09, and (d) <5 av =0.00, 
<5ip=-0.10, 5op = 0-10- The averaged hole-density, <5 av , is defined 
as, <5 av = (5ip + <5op)/2. Two thick (thin) lines correspond to 
the even (odd) symmetry band as regards the inversion. The 
doted line shows a half of Brillouin zone as a guide for eyes. The 
symbols, \n) (n=l-4), will be explained below. 
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Fig. 4. Energy scheme of four-layers composed of OPi, OP2, IPi 
and IP2- The energy of OPi and OP2 is separated from IPi and 
IP2 by W. Considering the inversion symmetry, these planes are 
classified into even- and odd wave function as OP C von, OP 0( jd, 
IPeven and IP dd- Including t±i, the degenerated states of IP OVC n 
and IP dd are separated. 



